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We present the results of a search for the radiative decay B — > Kn'j and find evidence for 
B + — > K + n / 'y decays at the 3.3 standard deviation level with a partial branching fraction of (3.6 ± 
1.2 ± 0.4) x 10 -6 , where the first error is statistical and the second systematic. This measurement 
is restricted to the region of combined Krf invariant mass less than 3.4GeV/c 2 . A 90% confidence 
level upper limit of 6.4 x 10~ 6 is obtained for the decay B° — > K°rj '7 in the same Kvj invariant 
mass region. These results are obtained from a 605fb _1 data sample containing 657 xl0 6 B-B pairs 
collected at the T(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. 

PACS numbers: 13.25.Hw, 13.20.He, 14.40.Gx, 14.40.Nd 



Radiative B meson decays proceed primarily through the flavour changing neutral current (FCNC) quark-level 
process b — + sj. FCNC processes are forbidden at tree-level within the Standard Model (SM), and hence b — > sj must 
proceed via radiative loop diagrams. As loop processes may include unknown heavy particles mediating the loop, any 
disparity between experimental measurement and SM prediction could be evidence of such new particles. 

The world average experimental branching fraction (BF) for the meson-level process B — > X s -f ((3.55 ±0.26) x 10~ 4 
Q) and the theoretical SM predictions ((3.15 ±0.23) x 10~ 4 0) are consistent. Measurements of individual exclusive 
B —> X s j modes, such as B — > Krj'j, provide consistency checks on the agreement between theory and experiment, 
and improve our understanding of the hadronization process in B — > X s "f and B — ► X s P~l~ . 

The analysis of the decay mode B — ► Krf^f uses 605 fb -1 of data collected at the T(45) resonance with the Belle 
detector. This decay was previously studied by BaBar who set upper limits (ULs) of B(B + — > K + rj'^j) < 4.2 x 10~ 6 
and B{B° K°rj'j) < 6.6 x 10" 6 at 90% confidence level (CL) from an analysis of 211 ftr 1 of data [|. 

A previous comparison of the modes B — > Kr\ and B — » Krj 1 showed a surprising suppression of the former 
with respect to the latter, due to destructive interference between two penguin amplitudes Q. The ability of QCD 
factorization techniques to correctly predict this BF hierarchy has been demonstrated, though the errors are large (fjj . 
A similar comparison of the observed B — ► Krj^f mode [(| 0] and the previous upper limits for B — ► Kr/'j shows the 
opposite trend for these BF's. The analogous QCD calculation has not yet been performed for these decay modes. 
Measurement of B — ► Kr/'^y is an important test of such a calculation. In addition, a time-dependent charge-parity 
(CP) asymmetry measurement will be possible if there are sufficient statistics for the decay B° — > K s rj'^ pf. Such 
mixing-induced CP asymmetries are suppressed within the SM, but some beyond-SM theories involving right-handed 
currents allow them to be large, even when the B° — > K s T)'~f BF agrees with SM predictions [Si. 

The Belle detector is designed to identify and measure particles from T(45) — ► BB decays [9|. It is located at the 
interaction point (IP) of the KEKB accelerator in Tsukuba, Japan, which collides electrons and positrons at energies 
of 8 GeV and 3.5 GeV, respectively [loj . Charged particle tracking and momentum measurements are provided by 
a silicon vertex detector (SVD) and a helium/ethane central drift chamber (CDC). Particle identification (PID) 
is performed using the CDC in conjunction with an array of aerogel Cherenkov counters (ACC) and time-of-flight 
scintillators (TOF). Electrons and photons are identified and their energy measured by an electromagnetic calorimeter 
(ECL) of CsI(Tl) crystals. These detectors are within a 1.5T magnetic field provided by a superconducting solenoid. 
A layered iron and resistive plate counter detector outside the solenoid gives Kl/^ separation. 

The primary signature of a B — > Kr/'j event is a high energy photon, which we reconstruct as an isolated shower 
in the ECL barrel (32° < < 129°), with shape consistent with a single photon hypothesis and no associated charged 
track. Only photons within the range 1.8 GeV < E* < 3.4 GeV are considered, where the asterisk denotes the e + e~ 
center-of-mass (CM) frame. All other photons in the analysis are required to have energies greater than 50MeV. 
Contamination from r] and 7r° mesons decaying to 77 is reduced using a likelihood technique based on probability 
density functions (PDFs) from Monte Carlo (MC) simulated data. If more than one photon in an event passes the 
selection, the one with the highest energy is taken. 

Candidate rf mesons are reconstructed in the r\' — > 7^7r + 7r _ and rf — > /9°7 modes. Here, r\ mesons are reconstructed 
as r\ — > 77 and r\ — > 7r + 7r~7r°, p° and Kg mesons are reconstructed as 7r + 7r _ pairs, and 7r° mesons as 77 pairs. 
Charged kaons and pions are separated using a PID selection with an efficiency for kaons (pions) of 85% (98%) and a 
misidentification probability of 12% (14%). PID information combined with a measurement of the energy deposited 
in the ECL is used to reduce electron contamination. 
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(c)Mbc in A_B signal region for B° — » (d)AE in Mb c signal region for B° — > ifSjj'7 

FIG. 1: Projections to the signal region from the fits to 605fb _1 of data. The solid lines are the combined background plus 
Krj'j functions, the dotted lines are the e + e~ — > qq functions, the dot-dashed the b — > c functions, the long dashed the Kr/'"/ 
functions and the short dashed the b — > u,d,s functions. 



K s candidates must pass a set of momentum-dependent selection criteria based on proximity to the IP, flight 
length, and the angle between the momentum vector and reconstructed vertex vector. Candidates with invariant 
masses lOMeVc 2 or more from the nominal K s mass of 497.7 MeVc 2 [lj| are rejected, where the allowed range is 
equivalent to a window of 4 standard deviations (a) of the natural mass width convolved with detector mass resolution. 
Neutral pion candidates must have momenta greater than lOOMeVc and 77 combined invariant masses in the range 
119MeVc 2 (2.5a) < M 77 < 152MeVc 2 (3a). 

rj — > 77 candidates are reconstructed using photons of energy greater than 100 MeV, and must have invariant 
masses in the range 490 MeVc 2 (4a) < M 77 < 590 MeVc 2 (3a) . They are also required to satisfy the helicity angle 
requirement | cos 6> hcl | < 0.9, where #hei is the angle between the momentum vectors of the 7/ and one of the decay 7's 
in the 77 rest frame. 77 — > 7r + 7r~7r° candidates must have invariant masses in the range 536 MeVc 2 (3a) < M^+^-^o < 
560 MeVc 2 (3a). The momenta of all 77 candidates are then corrected using a mass-constrained fit. p° candidates 
must have invariant masses within the range 550 MeVc 2 (3ct) < M w + W - < 950 MeVc 2 (3a) and also pass a helicity 
requirement similar to the 77. A fit constraining the 7r + 7r~ vectors to a common vertex must succeed. 

77' — * p°7 candidates must have invariant masses in the range 945 MeVc 2 (2a) < M p o 7 < 970 MeVc 2 (2a), and 
the photon must have energy greater than 200 MeV. ?/ — ► ?77r + 7r~ candidates must have invariant masses within the 
range 950 MeVc 2 (2a) < M V7r + v - < 965 MeVc 2 (2a). Both types of candidates are required to have momenta greater 
than 1.0 GeVc as measured in the CM frame. 

The invariant mass of the -^77' system (Mx-q 1 ) is required to be less than 3.4GeVc 2 . This requirement removes 
some background while retaining most of the expected Mk^ 1 spectrum. Two kinematic variables are defined: M^ c = 

(l/c 2 )^i? 2 oam — (p* Kr]l +p^) 2 c 2 and AE = E* B — E'bcam, where E* B is the energy of the candidate B meson in the CM 

frame, -Eboam is half the total CM energy (y/s/2), and p* Kv > and are the CM frame momenta of the Kr]' combination 
and the signal photon, respectively. In the calculation of Mb c , the momentum of the signal photon is rescaled to be 
p* = (l/c)(E hcam - E* Kn ,). B candidates must satisfy \AE\ < 0.3 GeV and 5.20 GeVc 2 < M hc < 5.29 GeVc 2 , and 
the signal region is defined as -0.1 GeV < AE < 0.07 GeV and 5.27 GeVc 2 < M hc < 5.29 GeVc 2 . 

Backgrounds to Krj'^j events are estimated using large MC samples (T^ |. The dominant background is from e + e~ — > 
QQ (q = u,d,s,c) continuum processes. To reduce this we form a Fisher discriminant [13\ from 16 modified Fox- 
Wolfram moments [T3| and the scalar sum of the event transverse momenta. A likelihood ratio (LR) is formed from 
the optimized Fisher discriminant, the cosine of the angle between the B meson flight direction and the positron beam 
axis, and the distance along the positron beam axis between the two B meson vertices, which are reconstructed from 
the SVD and CDC response to charged particles. The signal regions of the LR distributions are chosen by maximizing 
the figure of merit defined as A( s /(-»/A/ r g +A/" SB ), where Af s is the number of B — > Ktj'j MC events that lie above 
a certain LR value and Af SB is the corresponding number of sideband data events lying above the same value. The 
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TABLE I: The yields, efficiencies (e), daughter branching fraction products (FJ), measured branching fractions (B), signal 
significances including systematics (<S) and 90% CL ULs for the measured decays. 





Mode 


Yield (events) 


e 


n 


B(HT 6 ) 


S(a) 


UL(10~ 6 ) 


B+ 
B° 


-» K + rf-y 


°^- u -10.8 

c 1+5.0 
°- 1 -4.0 


0.024 
0.016 


0.571 
0.197 


3.6 ±1.2 ±0.4 

9 r + 2.4+0.4 
^• u -1.9-0.5 


H CO 
CO CO 


5.6 
6.4 



data sideband regions are defined as M bc < 5.26GeV/c 2 and either AE < -0.2 GeV or AE > 0.1 GeV. The BF 
central values found in the BaBar Kr/'j analysis [|[ are used to calculate J\f s , and the sideband data is scaled by 
the ratio of qq MC events in the full fitting region to those in the sideband region. To maximize discrimination, the 
optimization of th e fig ure of merit is performed in bins of flavour-tagging quality, as calculated by the Belle £?-flavour 
tagging algorithm [la ]. The LR requirements are 38% efficient for signal MC and remove 98% of background. 

A significant proportion of the MC events modelled as b — > c processes that pass the above selection criteria are 
found to include a D° meson. Any kaon candidate which forms an invariant mass within the range 1.84 GeV/ c 2 < 
^k-tt+ < 1.89 GeV/ c 2 when combined with any charged pion in the same event is removed from consideration, in 
order to suppress D° — ► K~ir + decays. B — > J/ijjK — ► {rj lr y)K events are suppressed by vetoing candidates with a 
combined 7/7 invariant mass within ±25McVc 2 of the nominal J/ip mass [llj . This reduces this background to a 
negligible level. 

On average, 1.24 candidates per signal MC event pass the selection criteria. A series of selection criteria are used 
to choose the best candidate, including the lowest B vertex % 2 , the lowest p° — > 7r + 7r~ or rf — > r\-K + -K~ vertex % 2 , the 
reconstruction with r\ candidate invariant mass closest to the nominal value [111 ], the highest E 1 from rj' — > p°7, and 
the lowest Kg vertex \ 2 ■ This technique selects the correct candidate in 76% of cases. 

Signal yields are extracted using extended unbinned maximum likelihood fits to AE and M^ c . All reconstructed 
charged final states are combined into one fitted distribution, and all reconstructed neutral final states are combined 
into another. PDF parameters are determined from MC distributions. The Krf^ distribution is modelled with a 
Crystal Ball line shape (CBLS) [l6| function for Mbc and CBLS plus Gaussian functions with common means and 
relative widths for AE. The b — > c background is modelled with an ARGUS [TtJ function for Mb c and a second order 
Chebyshev polynomial for AE. The b — > u,d,s background is modelled with a 2D Keys PDF [l8j]. The e + e~ — > qq 
distribution is modelled with an ARGUS function for M bc and a first order Chebyshev polynomial for AE. 

The means and widths of the CBLS functions describing the Krj'j distribution in AE and Mb c are calibrated using 
large control samples of B — > K*(892)^ data and MC. The PDFs of the Kr/'j, e + e~ qq, b c and b — ► u, d, s 
distributions are combined and used to fit the 605 fb -1 of accumulated data. The normalizations of the b — ► c and 
b — ► u, d, s components are fixed to values expected from MC studies; the Krj'j and e + e~ — > qq normalizations and the 
e+e~ — > qq function parameters are allowed to float, except for the ARGUS endpoint, which is fixed to 5.29 GeV/ c 2 . 
Figure [T] shows the results of the B + — > K +r q'^i and B° — > K^rf^ fits to data, where the Mb c plots show a projection 
in the AE signal region and the AE plots show a projection in the Mb c signal region. 

Table [I] shows the measured yields and signal significances for the fits to data. We find 33li[ 2 events for the fit to the 
charged modes and 5^4 events for the fit to the neutral modes. The signal significance is defined as yf— 2ki(Zo /L max ) , 
where L max and Lq are the values of the likelihood function when the signal yield is floated or fixed to zero, respectively. 
The systematic errors described below are included in the significances by convolving the likelihood functions with 
Gaussians of width defined by the magnitude of the errors. The signal significances including systematic errors are 
3.3(7 and 1.3tr for the charged modes and neutral modes, respectively. 

Any bias in the fitting process is determined using pscudoexperiments of MC scaled to the yields returned by the 
fits to data. The e + e~ — > qq and b — > c MC components are generated from the shape of the PDFs and combined 
with fully simulated MC for the Krj'^f and b — » u,d,s components. The analysis of 1000 pseudoexperiments finds no 
significant bias. 

The signal MC reconstruction efficiencies are calculated as the number of Krf n / events passing the selection criteria 
divided by the number generated. To calibrate for Mk v ' dependence, the efficiencies in 10 bins across Mk V ' are 
weighted according to the efficiency-corrected background-subtracted Mk v ' distributions in data. The means of these 
weighted efficiencies are listed in the e column of Table [TJ 

The BFs are calculated from the signal yields, calibrated efficiencies, daughter BF products, and the number 
of B mesons in the data sample. Equal production of charged and neutral B meson pairs is assumed. We find 
B(B+ -> X+7/7) = (3-6 ± 1.2 ± 0.4) x 10~ 6 and B(B° -> K°7]'-f) = (2 .5±lj±° A 5 ) x 10~ 6 , where the first errors 
are statistical and the second systematic. The 90% CL ULs are found to be B(B + — > K+rf-f) < 5.6 x 10~ 6 , and 
B(B Q — > K°rj'j) < 6.4 x 10 -6 . The ULs are calculated by integrating the likelihood function with systematic errors 
included in the physically allowed BF region. The UL is then defined as the BF below which 90% of the integrated 
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likelihood lies. The BFs and ULs are measured within the reduced phase-space region Mk V ' < 3.4GeV/c 2 . 

The systematic uncertainties for the charged (neutral) BF include errors on the following processes: photon detection 
(2.8% (2.8%)); 7T° reconstruction (0.8% (0.5%)); K% reconstruction (0.0% (4.5%)); r] reconstruction (3.4% (3.6%)); 
charged track detection (3.8% (5.0%)); K + /ir + differentiation (1.4% (1.5%)); and the calculated number of BB pairs 
in the data sample (1.4% (1.4%)). The statistical uncertainty on the MC efficiency after calibration is 1.7% (1.9%). 
The data/MC LR selection efficiency difference is estimated to be 3.7%(3.7%) using aB-> K* (892)7 control sample. 
The vetoes on M v ' 7 around the J/ip invariant mass and on M K -^+ around the D° mass are found to affect the 
reconstruction efficiency by 0.2% (0.4%)) and 0.5% (0.0%) respectively. The bias study results have an uncertainty of 
1.5% (3.5%). The correction of the signal reconstruction efficiency from the Mk-q' distributions has an uncertainty of 
1.2% (1.5%). The effect of Kgrfj events being reconstructed as K +r q'^ and entering the incorrect fit distribution, and 
vice versa, gives a — 6%(+6%) efficiency uncertainty. All fixed parameters in the fit to data are varied by ±lcr (the 
b — > c and b — > u,d,s yields by 3a), yielding uncertainties of (+6.5%, —6.7%) for B + — > K + r)'j, and (+11.7%, —16.6%) 
for B Q K^r]'j. 

The existing measurements of B — > X s ^f and B — > both rely heavily on the accuracy of the X s hadronization 

model performed by the JETSET [lj| program 20] [21] 22]. A large sample of inclusive B — > X s j MC was generated 
with the X s system hadronizcd by JETSET according to the Kagan-Neubert model [23[ and the mass of the b quark set 
to 4.75 GeV/c 2 . From this sample, the model predicts the BF of B + — > K + r]'^ events in the region Mkt) 1 < 3.4 GeV/c 2 
to be (1.8 ± 0.2) x 10~ 6 , while we measure (3.6 ± 1.2 ± 0.4) x 10~ 6 . In addition, this model predicts B(B -> %ito 
be (8.2 ± 0.9) x 10~ 6 using the same sample, which can be compared to the measured BF of (7.9 ± 0.9) x 10 -6 [ll| . 
With the current statistics, the BFs for both B + — > K + r)'^ and B — > if 777 are consistent with the Kagan-Neubert 
model. 

In conclusion, we report the first evidence of the decay B + — ► K + r)'^ with a partial BF of B(B + — > K + rj'^) = 
(3.6 ± 1.2 ± 0.4) x 10~ 6 and a significance of 3.3er in the region M Kr) , < 3.4GeV/c 2 . We also set a 90% CL UL of 
B(B" -> K°rf-/) < 6.4 x 10~ 6 in the same M Kr) , region. 
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